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Abstract Stellar halos may hold some of the best preserved fossils of the 
formation history of galaxies. They are a natural product of the merging pro- 
cesses that probably take place during the assembly of a galaxy and hence 
may well be the most ubiquitous component of galaxies, independently of 
their Hubble type. This review focuses on our current understanding of the 
spatial structure, the kinematics and chemistry of halo stars in the Milky 
Way. In recent years, we have experienced a change in paradigm thanks to 
the discovery of large amounts of substructure, especially in the outer halo. 
I discuss the implications of the currently available observational constraints 
and fold them into several possible formation scenarios. Unraveling the for- 
mation of the Galactic halo will be possible in the near future through a 
combination of large wide field photometric and spectroscopic surveys, and 
especially in the era of Gaia. 

Keywords Galaxy: halo • Galaxy: formation • Galaxy: evolution • Galaxy: 
kinematics and dynamics 



1 Introduction 



A word of caution is necessary before we start our journey through the stellar 
halo of the Galaxy. While working on this article, I could not help but wonder 
whether this was a good time to write a review on the stel l ar halo (see, for 
example, previous excellent reviews bvlGilmore et al.lll989l: iMaiews 5 1199.11: 
iBland-Hawthorn fc FreemanI I2OO0I : iFreeman fc Bland-Hawthornll20oi i. The 
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advances in the recent past have been enormous thanks to many (ongoing) 
large wide field photometric and spectroscopic surveys. In the past two years, 
almost every month a new satellite galaxy or a new stream was discovered. 
This in itself may be a justification for writing a review, but it also may lead 
to an article that is rapidly out of date. 

On the other hand, this of course, is only a sign of the field's health and 
its great promise for the young generations of scientists hoping to unravel the 
formation of the Galaxy. The reader is hereby warned that our understanding 
of the Galactic halo is still evolving significantly. In an attempt to avoid 
the (almost unavoidable) risk of becoming quickly outdated, the author has 
decided to review those properties of the stellar halo that (hopefully) provide 
the most direct clues to its evolutionary path. 

The stellar halo is arguably the component that contains the most useful 
information about the evolutionary history of the Galaxy. This is because the 
most metal-poor stars in the Galaxy and possibly some of the oldest ones are 
found here. It therefore provides us with a picture of the Milky Way in its 
early stages of evolution. Low-mass stars live for much longer than the present 
age of the Universe, and so retain in their atmospheres a record of the chemi- 
cal elements of the environment in which they were born. The very metal-poor 
halo stars are thus fossils, whose che mical abundance and motions contain 
information of their sites of origin. As lEggen. Lvnden-Bell k Sandaed (|1962| ) 
nicely put it: "... a study of these subsystems allows us partially to recon- 
struct the Galactic past because the time required for stars in the Galactic 
system to exchange their energies and momenta is very long compared with 
the age of the Galaxy. Hence knowledge of the present energy and momenta of 
individual objects tells us something of the initial dynamic conditions under 
which they were formed' . 

Despite its crucial importance, our knowledge of the stellar halo is frag- 
mentary. This is at least partly driven by the paucity of halo stars (near the 
Sun roughly only one star in a thousand belongs to the halo, as discussed 
in Section 12. 2p . This implies that generally a set of suitable (and preferably 
unbiased) selection criteria are used to find these precious fossils, as we shall 
see later in this review. 

The study of the formation of the Galaxy is particularly relevant at this 
point in time. There is a cosmological framework in place, the A cold dark 
matter (CDM) paradigm, which broadly desc ribes from first princ iples how 
structures in the Universe have evolved (e.g. ISpringel et al ] l200(il . and ref- 
erences therein). In this model, galaxies form hierarchically, through the 
amalgamation of smaller proto-systems. Despite the many successes of this 
theory, particularly on large scales, we do not yet understand in detail how 
galaxies form. Many of the shortcomings of this model have become evident 
through comparisons to the properties of o ur Galaxy and its nearest neigh- 
bors (e.g. the "missing satellites" problem. iKlvpin et al.lll999t iMoore et al.l 
Il999t ). Therefore, it becomes clear that studies of the Galaxy and its compo- 
nents are crucial for a proper and complete understanding of galaxy evolution 
in the Universe. 

This review describes the various characteristics of the stellar halo in the 
hope that these will aid us in our understanding of the formation of the 
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Galaxy. We start by discussing the spatial structure and global kinematics 
in Section [21 This approach is rather classical, since we shall discuss sepa- 
rately the presence of substructure in Section[3l In principle, substructure can 
strongly affect the determination of the density profile or of the shape of the 
stellar halo, e.g. by intr oducing overdensities which may not always be easy 
to id entify as such (e.g. iNewberg fc Yannyl 120061 : IXu et al ] 120071 : iBdT et al.l 
Given our current understanding of the Galactic halo, such substruc- 
tures may well dominate at large radii, but at small distances (i.e. the stellar 
halo near and interior to the Solar circle), the mixing timescales are suffi- 
ciently short that even if all of the stellar halo had been formed through a 
seq uence of dis sipationless mergers it would be completely smooth spatially 
fsee lHelmi fc'w hitc 199i). 

In Section 2] we describe our current knowledge of the metallicity distri- 
bution of halo stars, and particularly focus on what it reveals about star 
formation in the early Universe. We also discuss what chemical abundance 
patterns tell us about the process of formation of the Galaxy. Not only halo 
field stars are interesting, but also globular clusters and satellite galaxies give 
us insight into the characteristics and evolution of the Milky Way halo. These 
are discussed in Section [5l 

In Section [6] we discuss formation scenarios for the stellar halo taking 
into account most (if not all) observational constraints presented in earlier 
sections. Recent technological advances have allowed us to begin to probe 
halos in other galaxies. Their characteristics and a comparison to our own 
stellar halo are briefly mentioned in Section [71 Finally, we provide a brief 
summary and outlook in Section [8l 



2 Spatial structure and global kinematics 

2.1 Star counts: density profile and shape 

Star counts have been used in astronomy for over a century to derive the 
struc t ure of our Gal axy (and even of the Universe, e.g. lKaptevn fc van RhiinI 
I1920I : IShaplevI ilQlSh . They have led to the conclusion that our Galaxy has 
several physically distinct compon ents, namely the thin and thick disk, the 
stellar halo and the bulge/bar (e.g. lBahcall fc SoneiralllQSOl : iGilmore fc ReidI 
[l983i)- 

The structure of the stellar halo is intim ately linked to how the Galaxy 
formed ( Freeman fc Bland-Hawthornll200^ . Different scenarios predict dif- 



ferent shapes and correlations with properties such as age, metallicity, etc., 
as we shall discuss later. In the simplest case, the structure of the stellar 
halo could also give us insight into the structure of the dark matter halo. For 
example, if a significant fraction of the dark matter had been baryonic (com- 
posed by MACHOs. |Paczvnskilll986[ ). then the stellar and dark halos would 
presumably be indistinguishable. However, this scenario appears unlikely, as 
the microlensing surveys are unable to assign m ore than 8% of the matter 
to com pact dark objects ( Tisserand et al.l [20071 although see lAlcock et al.l 



(l200flf ) 'or a different while earlier result). 
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In modern cosmological models (namely, cold dark matter dominated) , it 
is now clear that the structure of the stellar and dark halo do not necessar- 
ily follow one another. The dark halo is much more extended, and probably 
follows an NFW profile (Navarro, Frcnk & Whit^ 19961). Its half-mass ra - 



dius hes around 150 kpc (jKlvpin et all 120021 : iBattagha et all l2005l HqoI), 
and therefore its struc ture is probably rel ated to the present-day nearby 
large scale tidal field ( Navarro et al.l [2004 ). On the other hand, the stel- 



lar halo is highly concentrated (e.g. the effective radius as traced by t he 
halo globular clusters is well within the Solar circle, iFrenk fc Wh"it3 (119821)) 



a property that likely reflects a very early formation epocri^l (|Moorell200ll: 
iHelmi. White fc ST)ringe]|[200l) . 

In the early years, star counts were based on photographic images, and 
were usu ally restricted to just a fe w fields. This was also the basis for the pi- 



oneering [B^cair^^SoneiS ( 1980l ) model of the Galaxy. These authors took 



the approach that our Galaxy is "just" like any other galaxy, and so it should 
be possible to reproduce the star counts using models that fit the light distri- 
bution of external galaxies. For example, they proposed the halo followed a 
"de Vaucouleurs" spherical profile, which is characteristic for elliptical galax- 
ies (although it is now known that not all elliptical galaxies are well-fit by 
this function; furthermore Morrison (1993) showed that a power-law p cx r" 
with n ~ —3.5 as proposed bv iHarrisI (|l976f) on the basis of halo globular 
clusters, provides a better fit for the stellar halo t han an R^^'^ law). 

Already in early studies ( Kinman et al.l [T966[ ) there was evidence that 



the stellar halo's shape varied with distance from the Galactic center. These 
studies favored a halo that was rounder in the outskirts (mino r-to-major axis 
ratio q ^1) and more flattened in the inner parts {q ^ 0.5; ISchmidtlll956l : 
lHartwicklll987t iPreston et al.lll99ll ) . Typically RR Lyrae and blue horizontal 
branch stars (BHB) were used as tracers. Even today these stars are com- 
monly used to isolate samples of true halo stars. Such stars are only found in 
old populations, and they are relatively bright so they can be observed with 
mod est telescopes at dis tances as large as 100 kpc from the Galactic center 
(e.g. IClewlev et al. I I2005D . 

The density proflle of the stellar halo is often parametrized in a principal 
axis cartesian coordinate system as 

2 p \«/2 

p[x,y,z)^po (1) 

where n is the power-law exponent (and because of the finite extent of the 
Galaxy, it is negative), and q and p are the minor- and intermediate-to- 
major axis ratios (in the axisymmetric case, p = 1). The scale radius a is 
often neglected because o f the rather concen t rated nature of t he stellar halo 
(although see for example IPhleps et al.ll2005t iBica et al.ll2006[ ). Finally, po is 
the stellar halo density at a given radius tq, which is generally taken to be the 
solar radius for obvious reasons (and referred to as the local normalization). 

^ This refers to the inner stellar halo; as we shall show in Section |3] it is clear 
that the outer stellar halo is still forming. 
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As we shall see below, this parametrization is likely to be an oversimplification 
because the shape as well as the exponent seem to vary with distance from 
the Galactic center. 



2.1.1 Modern pencil beam surveys 



in et al. I (|2onnh have used a large combined set of deep star counts at 



high and intermediate galactic latitudes, and foun d a power-law index n 
-2.44 and flattening q = 0.76. ISiegel et al] (|2002D observed seven Kapteyn 
selected areas (14.9 deg^) (hence ensuring uniform photometry) and preferred 
n ~ —2.75 and q — 0.6. There are many delicate points in these analyses. 
For example, the counts also include stars from the thin and thick disks, 
which implies that the fits should also allow for the (self-consistent) modeling 
of these components (see top panels in Figure [T]) . Furthermore there is a 
degeneracy in the modeling of the stellar halo in the sense that a less flattened 
halo requires a steeper density proflle (bottom panel of Figured]). However, it 
is worth bearing in mind that in both studies the fits obtained are not perfect, 
and systematic differences between models and data ar e apparent, especially 
at the faint end (see also Figure 6 in lRobin et al.l (2003))- In particular, there 



are more stars towards the anticenter than predicted, and there are less stars 
observed towards the center. This is similar to what orrison et al. 1 (120001) 



find using turn-off halo stars (i.e. those with (B — V) ^ 0.38), who prefer n = 
—3 and q = 0.6. These authors note that a two-component halo is unlikely 
to solve these problems, because it is incapable of enhancing the number 
counts towards the anticenter (a rounder halo would in fact tend to lower the 
counts). Some proposed solutions are that the thick disk scale height varies 
with radius, or that the outer halo is lumpy. Either may well be possible. 
For example, fields at Galactic latitudes |6| < 25deg are now known to 
contain sign ificant substructure towards the anticenter ( Newberg et al. 200^ 
llbata et al.il2003l: lGrillmaiij|2006ah. M ore recently CADIS (Calar Alto Deep 



Imaging Survey; iPhleps et al.l (|2005h ) mapped 1/30 deg^ in the magnitude 



range 15.5 < i? < 23, and found n = —2.5 if g = 0.6. 



2.1.2 Panoramic surveys 

The advent of large galaxy surveys has been extraordinarily beneficial for 
Galactic astronomy. Although initially conceived to be a by product, multi- 
color deep wide field surveys have yielded completely new views of the Galac- 
tic halo. Some of the most spectacular results come from the Sloan Digital 
Sky Survey (SDS©. The first studies o n the Early Data R elease (EDR; 100 
deg^ down to r ^ 2 0) were carried out bv llvezic et al ] (|200nt) using RR L yrae; 
lYanny et al.l (j2000l ) using A-stars (BHBs and blue stragglers) and I Chen et al.l 
1"^ 



(j200l[ ) with stars of halo turn-off color. In agre ement with previou s work us- 
ing the same technique (see previous section) , IChen et alT ()200ll ) measured 
n = —2.3/2.5 a.nd q = 0.6/0 . 5 dep ending on the local normalization of the 
halo. However. lYannv et al.l ( 20001 ) estimated n = —3.2 ± 0.3 and q — 0.65 



, http://www.sdss. org, 
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Fig. 1 Top panels: Color distributions of two fields used bv lRobin et al.l (|2000l ) to 
constrain the sliape of the stellar halo. The dots correspond to the true number 
of observed stars in each color bin, while the heavy solid lines are the number of 
stars predicted by the model assuming a power law index of —2.44, a flattening 
of 0.76 and a local density of 0.75 xpo for the stellar halo. Thin lines show the 
contribution of the halo alone. The red peak corresponds to nearby M dwarfs, 
while the blue peak indicates the location of the halo and thick disk turn-off points. 
Bottom panel: Iso-likelihood contours at 1, 2 and 3 sigmas in the plane (power law 
index, flattening) of the spheroid population, coded by the local density, going from 
1 (solid), 0.75 (dotted); 0.50 (dashed) and 1.25 (dash-dotted) xpo- This diagram 
evidences the strong covarianc e betw een the shape and power-law index of the 
stellar halo. [From Robin et al] l|2000l) . Courtesy of Annie Robin. Reproduced with 
permission of A&A.] 



(after maskin g out the overden sities due to the Sagittarius streams), in accor- 
dance with llvezic et al.| ()200Cl[ ) who found n ~ —3. Note that the latter two 
works employ tracers that probe considerably larger distances (particularly 
in comparison to the halo turn -off stars, which are intrins ically much fainter: 
Mvt:o ^ 4 versus Mvhb ~ 0.6: lBinnev fc Merrifieldlll998t) . 

Using RR Lyrae stars, the QUEST survey (|Vivas et al.ll2004D find n = 
—3.1 ± 0.1 (and that s teeper profiles produ ce worse fits), and q varying with 
radius as proposed by iPreston et al.l (|l99lh . 
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The SDSS has no w imaKed 8000 deg^ of the sky yiel ding; a truly panoramic 
view of the Galaxy ( Adelman-McCarthv et anl2007t l. However, instead of 



simplifying the picture and providing definitive answers, it has shown that the 
halo is far more complex than envisioned some ten years ago. For example, 
iNewberg fc Yannvl (|2006h find that there are significant asymmetries with 



respect to the line connecting the galactic center to the anticenter in the 
number counts of halo turnoff star s. They postulat e this could be due to a 
triaxial halo. In concert with this, IXu et al. 1 (|2006D find an excess of stars 
to wards \l,b] = (270° 60°) compared to [l,b] = (90°, 60°). 

I Juric et "aP ( 2008| l have derived photometric parallaxe^ for 48 million 



stars in the SDSS catalog and attempted full modeling of the various Galactic 
components. These authors find that it is possible to measure the properties 
of the thin and thick disk relatively robustly (provided large portions of 
the celestial sphere have been observed). However they acknowledge that 
determining the properties of the halo is more difficult, partly because of 
the presence of significant clumps and overdensities, as shown in Figure [2l 
The halo power law index n is poorly constrained by the data, but an oblate 
halo with q — 0.5 seems to be favored, alth ough no good fits are obtained at 
the faint end. Amongst various structures. lJuric et al.l ( 2008( ) have detected 



a huge overdensity towards the constellation of Virgo at I = 270° and b 
60° (the same region that cannot be reproduced by smooth axisymmetric 
models). 

The SDSS is focused on the northern galactic cap. This implies that the 
nature of the overdensities (global versus local) may be difficult to establish 
without good coverage of the southern hemisphere. In particular, if the halo 
were triaxial, one would still ex pect (in the mos t reasonable case) symmetry 
with respect to the z = plane. iXu et al.l (2007) carried out this exercise and 
compared star counts towards the North Galactic Pole (NGP) from SDSS and 
towar ds the South Galact ic Pole (SGP) using SuperCOSMO^ photographic 
data (jHamblv et al.ll200lh . They found that there is an asymmetry between 



the north and southern cap of approximately 17% in the _B-band. Therefore, 
they concluded that the Virgo overdensity is likely a foreign component of 
the stellar halo (see also Section [3]). Their modeling, after subtracting the 
overdensities, favors an axisymmetric halo with parameters n — —2.8 and 
q = 0.7. 

This last example shows particularly well the need for large sky cover- 
age and deep photometry to establish the structure of the stellar halo of the 
Galaxy. This is the only way to break the known degeneracies affecting the 
modeling of number counts from restricted pencil beam surveys. Neverthe- 
less, given the problems in reproducing star counts at the faint end and the 
discovery of large overdensities on the sky, it is essential that the presence of 
substructure is taken into account in the models. 

At this point in time, it transpires that the inner halo may well be axisym- 
metric and flattened {q ~ 0.6), w hile the outer halo (t raced by e.g. RR Lyrae 
stars and very deep counts as in iRobin et al.l ( 20071 )) may be more spheri- 



^ assuming a "mean" photometric parallax relation, as function of color index, 
for all metallicities. 

,http://www- wfau.roe.ac.uk/sss/ 1 
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Fig. 2 Left panels: Stellar number density as a function of {R, z) for different 
r — i bins using SDSS photometric parallaxes (the distance scale greatly varies 
from panel to panel). Each white dotted rectangle denotes the bounding box of 
the region containing the data on the following panel. Note the large overdensities 
particularly in the top panels. Right panels: The stellar number density for the 
same color bin as in the top left panel (0.10 < r — i < 0.15, hence preferentially 
selecting stars near the halo turn-off point), shown here for cross-sections parallel 
to the Galactic plane, as a function of the distance from the plane (from 16 kpc 
at the top left to 6 kpc in the bottom right). The origin marks the Sun's position, 
while the circles aid in visualizing departures from axial symmetry of the Galaxy. 
Note the strong asymmetry with re s pect t o the y = G line, which is due to the 
Virgo overdensity. [From Jur ic et al.l l|2008l ). Courtesy of Mario Juric. Reproduced 
with permission of the AAS.] 



cal (and steeper). It would be interesting to explore whether this is due to 
different formation paths for the inner and outer halo (e.g. dissipative ver- 
sus non-dissipative formation, see Section [6TT|) . The gravitational pull of the 
stellar disk (and the adiabatic contra ction of the halo) may also be (partly) 
responsible for this dichotomy, see e.g. lBinnev fc Mavl ( 1986f ): [Chiba fc BeersI 
([2001). In principle, it is also possible that the dual shape of the stellar halo 
reflects a change in the number of progenitors and their orbital characteristics 
with distance from the Galactic center. 



2.2 Normalization 



The local normalization (density of the stellar halo at the position of the Sun) 
has been determined both locally as well as from deep pencil beam surveys. 
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Im orriso isolated local halo giants from an objective-prism survey of 

metal weak stars, and derived a halo number density expressed in terms of 
halo-to-disk ratio of 1:1200 for K giants (My < 0.5), and 1:850 for dwarf 
stars with m agnitudes 5 < My < 8. The latter is comparable to the estimate 
obtained bv lFuchs fc Jahreil ( 19981 ) using nearb y M dwarfs (withi n 25 pc), 
who measur ed pn = 1-5 x 10""^ M f^/pc^ (see also lDigbv et alJl2003[ ). On the 
other hand, iFlvnn et al.l ( 1996t ) used very deep HST counts of M dwarfs 
and obtained a factor four smaller number of red dwarfs, in agreement with 
iGould et all (|l998. ) who find po 6.4 x 10"^ Mg/pc^ using subdwarfs in 53 
fields observed with the Wide Field Planetary Camera on the Hubble Space 
Telescope (HST). Attempts to measure the normalization using BHB stars 
are more tricky because of the variation of the properties of the horizontal 
branch with age and metallicity. This implies that the transformation from 
number of BHBs to Mq is non-trivial. Note that the measurement of po 
relies on an accurate determination of the luminosity function (and hence 
on the slope of the initial mass function), but possibly the largest source of 
discrepancy is the varying shape of the halo with radius. For example, the 
deep counts using HST explore larger distances, and hence presumably also 
the rounder part of the stellar halo, and therefore necessarily favor a rather 
low local normalization (jFuchs &: Jahreifi,il998i ). 



2.3 Kinematics 

Obtaining space velocities for halo stars is no doubt more laborious than map- 
ping their sky distribution. Although radial velocities are relatively straight- 
forward to obtain for a small number of objects even for rather faint stars 
(e.g. y ~ 19 take less than one hour-long integrations with 4m-8m class 
telescopes) , measuring their tangential velocities is only possible if both the 
proper motion and distance are known. While estimates of the latter may be 
obtained through photometric or spectroscopic parallaxes, proper motions 
from the ground require a large time basel ine (of the order of 100 years) to 
produce miUiarcsecond accuracy (see, e.g. iBeers et al. 2000l and references 



therein). Therefore, and for many other reasons ([Lindegren 120041) . astromet- 
ric space missions with wide field angle capabilities are an absolute necessity. 

Presently the catalogs of halo stars with relatively accurate 3D kine- 
matics are somewhat small, and typically sample a few kpc around the 
Sun. They have been constructed using proper motions mainly from Hip- 
parcos, and radial velocities from other surveys. Velocity errors are in the 
ra nge 20 - 30 km s~ ^ . Currently the largest such catalog available is that 



bv iBeers et al. 1 (I2000D . with approximately one thousand nearby stars with 



[Fe/H] < -0.6 dex. 
2.3.1 Local halo 

Local samples of halo stars show a small amount of prograde rotation V^, ^ 
30-50 km s-^ (ICarnev et al.lll996l: iLavden et al1ll996l : IChiba fc Bee"^l2000l: 



Im orrison et al.ll 



2008ah The velocity ellipsoid is roughly aligned with a cylin- 



drical coordinate system, having dispersions of (crj^, (t^, Uz) — (141±11, 106± 
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9, 94 ±8) km (|Chiba fc Beerd[200nh . The stellar halo is therefore sup- 
ported by random motions, and its flattened shape is con sistent with the de- 

gree o f velocity anisotropy observed (see Section [2.3.3[ also lSommer-Larsen fc ChristensenI 
EMi). 

The velocity distribution of halo stars near the Sun appears to be fairly 
smooth, with only a small degree of substructure (see Section [3] for more de- 
tails). There is no strong correlation of the rotational velocity with metal- 
licity. However, at large distances ab ove the Galactic plane ( ~ 3 — 5 kpc), 
the ha lo become s slightly retro g rade ( Maiewski et all ( 19961 ): ICarnev et al.1 
1996D. see also iKinman et al.l (l2007f)). a nd its mean metallicity is lower 



Carollo et al.ll2007HMorrison et al.ll2008al ) 



2.3.2 Global properties 

Figure [3] show s the line of sight velocities of a sample of 240 objects in the 
Galactic halo (|Battagha et al.ll2005L [2006h as function of distance from the 



Galactic center. This sample includes K giant stars from the Spaghetti sur- 
vey, field horizontal branch stars, globular clusters and dwarf galaxies. On the 
right-hand side panel, the radial velocity dispersion is plotted. This quantity 
shows an almost constant value of 120 km out to 30 kpc and then continu- 
ously declines down to 50 km s~^ at about 120 kpc fsee lSommer-Larsen et al] 
(|1997D for an earlier similar conclusion from a smaller and less distant sam- 
ple) . This fall-off puts important constraints on the density profile and total 
mass of the dark matter halo of the Milky Way. In particular, under the 
assumption of a power-law density for the stellar halo and a constant veloc- 
ity anisotropy, an isothermal profile ca n be ruled out, wh ile both a dark 



halo following a truncated flat model ([Wilkinson fc Evans 1999) of mass 
5lij X 10" Mq and an NFW profile of mass 9.4lJ;^ x 10" Mfr, and con- 
centration c = 18 are consistent with the data (see also iMoore et al.l 120061 : 
lAbadi et al.ll2006D . 

The stellar halo is not completely smooth kine matically (see Section [3l. 
Evide nce of halo moving groups was presented by iRatnatunga fc FreemanI 



(I1985D towards the South Galactic Pole using radial velocities for K giants. 



bv 'M ajewski et all ( 199($ ) towards t he North Galactic P o le bas e d on proper 



motio ns, a nd more recently by e.g. IClewlev fc KinmanI ( 2006f ). ISirko et al] 



(I2004D and I Vivas etaH (|2004D using horizontal branch stars' radial velocities. 



We discuss these substructures in more detail in Section [31 



2.3.3 A small digression on the relation between the shape of the stellar halo 
and its velocity ellipsoid in the Solar neighborhood 

As described above, the stellar halo docs not show significant rotation. This 
implies that its flattened shape has to be attributed to anisotropy in the 
velocity ellipsoid [^ 



^ A similar argument i s used to explain the shape of elliptical galaxies, see 
IBinnev fc Tremaind ([19871 1. 



The stellar halo of the Galaxy 



11 




Fig. 3 Radial velocities in the Galactic Standard of Rest Vgsr as function of dis- 
tance from the Galactic center (left); and Galactocentric radial velocity dispersion 
of the Milky Way stellar halo (right). The curves show the best-fit models obtained 
assuming a truncated flat dark matter halo of mass S^i'y x lO^'^ Mq (solid curve) 
and an NFW profile of mass 9.4tJ * X 10" Mo and c= 18 (dotted curve), with 
constant anisotropy. The strong decline suggests that for an NFW profile to pro- 
vide a better fit eilSher the anisotropy must change with radius, or the stellar halo 
must have an edge. [From iBattaglia et al.l (|2005l . [20061 ) . Courtesy of Giuseppina 
Battaglia. Reproduced with permission of Wiley-Blackwell Publishing.] 



We will now derive the expected relation between the flattening of the 
stellar halo and the degree of anisotropy using the virial theorem. In steady 
state, the virial theorem in tensor form may be expressed as 

2Tjfe + n,k + Wjk = (2) 

where 

with /3(x) the density profile of the system, Vi the mean velocity along the 
i-axis, (t|;. its velocity tensor, and ^(x) the gravitational potential of the 
dominant mass component. Note that in the case of the stellar halo, the 
density and the gravitational potential are not related through Poisson's 
equation; the stellar halo is a tracer population embedded in an underlying 
potential whose main contributor is the dark matter halo. 

If we assume that the dark halo is spherically symmetric and that the 
stellar halo is axisymmetric, then Wxx = Wyy and Wjk = for j ^ k. 
Similar arguments can be used to show that all equations are automatically 
satisfied for j ^ k, while the only two independent equations that remain are 



= 2 y d^xpVjVk 

f ^3 9$ 
= / d-^xpzh^ — 

J OXj 



+ + Wxx^O 
2T,, + 77,, + W,, = 0. 
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Because the streaming (mean) motion of the steUar halo is neghgible 
compared to the random motions, we may set T^x ~ T^z — 0. In this case, 



Wx 



Wz, 



(3) 



which states that shape of the mass-weighted velocity dispersion tensor de- 
pends on the shape of the (moments of the) mass-weighted gravitational 
potential. Equation ([3]) is equivalent to 



/ (Px p{R, z) (tI 
J (Px p{R, z) a1 



Jd^xp{R,z)-- 

z^ df^ ' 

f (Px p(R, z) — 

r dr 



(4) 



Note that ax 7^ cr^ if the density distribution of the tracer is axisymmetric. 
For simplicity we will assume 

1. The shape of the velocity ellipsoid does not depend on location 

(^l/al = constant. (5) 

2. A logarithmic potential for the dark halo 

<P{r) ^v^Jn{r^ +d'^). (6) 

3. The density profile of the stellar halo follows a power-law, is axisymmetric 
and has flattening q. Using Eq. [T]we obtain 



p{R, z) = po{m/m,or, nP = R^ + z^/q^ 



(7) 



Under these conditions, the left-handside of Eq. ([4]) reduces to a constant. 
The goal now is to evaluate the right-handside expression and compare it to 
the numerical value derived observationally for this constant. 

Let us consider the case in which n — —A. This value is somewhat steeper 
than that discussed in Section 12.11 but it has the advantage of allowing the 
integrals to be solved analytically, and hence leads to a simpler expression. 
In this case Eq. (|4]) becomes 



1 



4(7^(1 — (72)arctg 



q^/l - g2 + 2(1 - 2(7^)arctg 



1 + q 



(8) 

The observed flattening of g ~ 0.6 therefore translates into a predicted 
al/al = l/2{(jl + al)/(jl = l/2{ a\ + al)/al = 1.73. The measured value of 
l/2(cr| -h (Tl)/al = 1.76 ± 0.37 (|Chiba &: Beersll2000[ ) is interestingly close. 
This presumably also means that the effective Galactic pote ntial near the 
Sun is quite close to that of an is othermal sphere (see also I White! Il985t 
ISommer-Larsen fc Christensenlll989[ ). 
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Fig. 4 A "simulated" stellar halo built up from accreted satellites. The color scale 
indicates surface brightness. [Sanjib Sharma, Kathryn Johnston and James Bullock 
are acknowledged for this figure.] 



3 Substructure 



One of the most fundamental ingredients of the ylCDM cosmological model 
is that galaxies grow via mergers. Therefore, perhaps the most direct way of 
testing this paradigm is by quantifying the amount of mergers that galaxies 
have experienced over their lifetime. This implies finding the traces of those 
merger events. 

Mergers are expected to have left large amounts of debris in the present- 
day components of galaxies. In particul ar, a s tellar halo is easily built up via 



the superposition of disrupted satellites ((Johnst on et al. 1996; Helmi & WhitJ 
119991: iBullock et all l200ll : iHarding et all 1200 1[) . This is clearly exemphfied 



in Figure [H adapted from the simulations and semi-analytic modeling by 
IBullock &: Johnston! (I2005D . 



The prediction that debris in the halo should be ubiquitous in hierarchi- 
cal cosmologies, together with the discoveries of a disrupting dwarf galaxy in 
the halo by Ibata et al. ( 1 994) as wel l as of nearby remains from an ancient 
accretion event bv lHelmi et al.l (|l999l ). have boosted the search for substruc- 
ture in this component of our Galaxy in recent years. Such searches have 
benefited enormously from the advent of large photometric surveys, and the 
SDSS in particular, as we shall see below, have all led to a shift in paradigm 
in the field over the past 10 years. 
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Fig. 5 Simulation of the disruption of a satellite galaxy orbiting in the inner 
regions of the halo of the Milky Way. The top panels show the xz projection of the 
particles' positions, which evidence significant evolution over short timescales. The 
velocities of the particles located in a spherical volume of ~ 2 kpc around the Sun 
at t~ 8 Gyr are strongly clu s tered , as shown in the bottom panels. Based on the 
simulation from I Helmi et ahl (|1999l ). 



From a theoretical point of view, one expects substructures in the spatial 
domain, as well as in velocity space (or in phase-space). One can distinguish 
essentially two regimes in which substructures will be most easily found in 
either of these domains: 



t ^ torb Short after infall, or for satellites orbiting the outer halo (dy- 
namical timescales are comparable to the time since infall). In this case, 
streams are very coherent in space. Hence, a simple way to reveal these 
accretion events is t o map of the positions of halo stars on the sky 
(| Johnston et al.lll996f ). 

t » torb Long after infall, or in the inner halo (r < 15 — 20 kpc). The 
debris is spatially well-mixed, as shown in Figure[5l Because of the conser- 
vation of phase-space density, the stars in a stream will have very similar 
velocities (see bottom panels in Figure O. In fact, streams get colder as 
time goes by. Therefore the characteristic signature of an accreted galaxy, 
even many Gyr after infall, is the presence of very cold structures in ve- 
locity space (|Helmi fc Whitdll999l) . 
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3.1 In the spatial domain 



Because of the relative ease with which young substructures may be uncov- 
ered in the spatial domain, the largest fraction of known streams are those 
that are spatially coherent. 

These discoveries did not really happen systematically until the advent 
of deep photometric wide-field surveys. The Sagittarius ( Sgr) dwarf galaxy - 
was a true beacon, bec ause it was clearly tidally disturbed (jlbata et al.l 19941 : 
iGomez-Flechoso et al.|[l999,) Shor tly after its discove ry, a southern ex t ension 
was found by Mateo et al.l ( 1998h . Somewhat later, iMajewski et al.l ( 19991 ) 



reported debris at 40° further south from the main bo dy, which could be 
tentatively fit by models of a disrupting dwarf ((Johnston et al.|[l999t ). A 
very large excess of RR Lyrae and BHB star s at high Galactic latitudes was 
subse quently discovered in the SDSS survey ( Ivezic et al.ll2000l : lYannv et al.l 
I2OOO) . This overdensity was roughly along the north-south elongation of the 
Sgr dSph, but now at 60° away from it, and at a significantly larger distance 
(50 kpc instead of 15 kpc from the Galactic center). The association of the 
overdensity to the Sgr dwarf galaxy northe rn/leading strean i was nearly ob- 
vious (see e.g.l Gomez-Flechoso et al ] |1999[ ). as the models of lHelmi fc Whitd 
(I2001D showecS 

After this, the Spaghetti Survey also found an overdensity at t he same lo- 



catio n using red giant b ranch st a rs wit h measured radial velocities (jPohm-Palmer et al.l 
|2001[ ). iMartfnez -Delg ado et al.l (|200lh using deep color-magnitude diagrams 
det ected the presence of a main sequence and turn-off at similar locations, 
and llbata et al. I (|200lh using carbon stars were able to trace the stream over 
an extended portion of the sky. 

A spectacular view of the Sgr streams was given by the Two Micron 
All Sky Survey (2MAS^|). iMaiewski etal\ (200a ) using M giants (selected 
from the 2MASS database on the basis of their colors) were able to map the 
streams 360° around the sky as shown in Figure El The vast extent of the Sgr 
dwarf streams has been used to constrain the shape of the Galactic gravita- 
tional potential, albeit with conflicting results (|Helmill2004bl : I Johnston et al.l 
120051: iFellhauer et al.l[2006h . 

In the years to follow, many o f the structures discovered in the halo 
would be Unked to the Sgr d ebris (|Newberg et all 120021 : IVivas et al.ll200ll : 
iMartinez-Delga do et al. 20041). An exception to this is the Monoceros ring 
(|Yannv et aTW ooa: Ib ata et al.| [2"003). a low-latitude stream towards the an- 
ticenter of the Galaxy, that spans about 100° in longitude at a nearly constant 
distance. It is somewhat of a semantic distinction to state that this structure 
is part of the stellar halo, or of the Galactic thin or thick disks. As yet, the 
nature of the feature an d of the Canis Major overdensity in a similar location 
at low Galactic latitude ( Martin et al.|[20M ). are highly debated. Possible in- 



terpre t ations are that it is debris from an accreted s atellite (e.g.lHelmi et al 



2003bl : iMartin et al\ l2004l : iPefiarrubia et all 120051 : iMartmez-Dekado et al 
2005), an overdensity associated to the Galactic warp (|Momanv et al. 1 12004 



® This was posted as a preprint on astro-ph in March 2000; just as the SDSS 
discoveries were being made. 
^ ,http://pegasus.phast.umass.edu7] 
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Fig. 6 2MASS revealing the streams from Sgr. Smoothed maps of the sky for point 
sources selected according to 11 < Ks < 12 and 1.00 < J — Ks < 1.05 (top), and 
12 < Ks < 13 and 1.05 < J — Ks < 1.15 (bottom). Tw o cycles around t he sk y 
are shown to demonstrate the continuity of features. [From Maicws ki et al.l (|2003l ). 
Courtesy of Steve Majewski. Reproduced with permission of the AAS.] 




Fig. 7 The field of streams from SP SS featuring Sgr, severa l new satellites and 
the low- latitude Monoceros ring. [From lBelokurov et al.l H2006bl ). Courtesy of Vasily 
Belokurov. Reproduced with permission of Wiley-Blackwell Publishing.] 



|2006[ ) , or even a proj ection effect caused b y looking al ong the ne a rby No rma- 
Cygnus spiral arm ( Moitinho et al.l l2006t ) . Recently, 'Grillmair ('2006a) has 
mapped the structure of the Monoceros overdensity using SDSS, and showed 
that it consists, at least in part, of a set of very narrow low-latitude sub- 
structurc s (see F igure [71), which can easily be explained by the models of 
[Penar ru bia et al.l ^2005). It is worthwhile mentioning here that such a satel- 
lite may well perturb the Galactic disk, and induce departures from axial 
symmetry or overdensities such as e.g. Canis Major. Furthermore, the Sgr 
dwarf leading tail also appears to be crossin g the Galactic plane at roughly 
the same location and distance as the ring f' Newberg et al.ll2007t) . and may 
well be responsible for some of the substructures observed. 
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At the time of writing of this review, other su bstantial substru ctures 
have been uncovered, such as the Virgo over densitv (IJuric et al.ll20"08h . pos- 
sib ly linked to the Vi rgo Stellar Structure (|Vivas et alJ ( 2001 ); confirmed 
bv iDuffau et alJ (l2006l) using radial velocities), the Triang ulum- Andromeda 
(iRocha-Pinto et al.l 2004h and the Hercules- Aquila clouds (jBelokurov et alJ 
l2007bD . The nature of these substructures is unclear, particularly because of 
the large area on the sky that they cover. It is evident that detailed kine- 
matic studies of large samples of stars with good distances, as well as chemical 
abundances, and full coverage of the structures, are necessary to understand 
whether these are truly related to tidal debris. In this case, one would expect 
to see coherent motions, as indeed appears to be the case for the Hercules- 
Aquila cloud. 

Narrow streams have been detected around glob ular clusters (iGrillmair et all 

Il995[ ). the mos t spectacular examples being Pal 5 (lOdenkirchen et al.l 2001 ) 
and NGC5466 (|Grillmair fc .Tohnsonl200fil : lBelokurov et al.l2006ah . This shows 
that gl obular clusters probably also have contributed to the halo field popu- 
lation ( Baumgardt et al. 2008t) . Another example of a nar row tidal tail is the 
Orphan-Stream (lGrillmair|l2Q06b|; iBelokurov etaHl2007ct) . whose progenitor 
may well have b een a d warf spheroidal-like object accreted by the Galaxy 
several Gyr ago (j Sales e t al. 2008). 

It is very important to understand what kind of biases affect our ability 
to detect substructure, particularly if our goal is to determine the mass func- 
tion of the objects accreted and the timing of these events. For example, it 
is obviously easier to detect high surface brightness features like those orig- 
inating in massive recently accreted satellites. For instance, the Sagittarius 
stream re aches a surface brightness of E y ~ 28 mag/arcsec^ at I ^ 350° and 
b 50° ('Martfnez-De lgado et al.l|200ll ). where it was first detected, while 
the other overdensities discovered so far typica lly have surface brightn esses 
Sv ^ 32 — 33 mag/ arcsec^ (see Figure 8 in IBelokurov et al ] l2007aL and 
private communication). 

The large majority of the above-mentioned substructures are located in 
the outer halo, at distances >15 kpc from the Galacti c center. The inner halo, 
on the other hand, appears to be relatively smooth ( Lemon et al.|[200^ . 



3.2 In the velocity domain 

As briefly discussed in Sec. 12.31 kinematic information is more expensive to 
obtain, and hence the velocity domain has been less explored despite its 
significantly larger suitability for searches of substructure. 

Two new major spectroscopic surveys are likely to change this in the 
nearby future. Th e Sloan Extension for Galactic Understanding and Explo- 
ration (s egue!!!, iBeers et al.l l2006l) and the RAdial Velocity Experiment 



(RAVe!3, ISteinmetz et all l20(i(Tl)~ are both currently collecting spectra for 



very considerable samples of stars. The final catalogs will be highly com- 
plementary. SEGUE will observe ~ 240, 000 stars in the range 15 < V < 21, 



9 



http:/ /www.sdss.org/segue/a boutseg ue.htm 1 1 
http://www.ra ve-s u rvey.aip.de/rave/| 
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while RAVE aims at 10^ stars with 9 < / < 12. The average velocity er- 
rors that these surveys can achieve are of the order of 10 km s~^ and 1 
km s~^ respectively. The advent of such large spectroscopic surveys is clearly 
a very interesting development since the kinematics of halo stars can be used 
to find and constrain the nature of the o yerdensities, as we ll as to measure 
the mass distribution in our Galaxy (e.g. ISmith et aLlfeoOTl ). 

The catalogs of halo stars with 3D kinematics used so far in searches 
of substruct ure are significantly smaller and less accurate (see for example, 



the work of 



Chiba fc YoshiH [1998 : Helmi et al.lll999t IChiba fc Beera .200d 



iKeplev et al.ll2007^ . Thanks to the developments mentioned above, however 



we may expect soon a ten-fold increase in the number of nearby halo stars 
with full kinematics. For example it will be possible to combine the RAVE 
radial velo cities with the relativ ely accurate proper motions from the Tycho- 
2 catalog. ( CaroUo et al.l (|2007t ) have already demonstrated the potential of 



combining the SDSS data (with radial velocities and stellar parameters) with 
proper motions from the USNO-B catalog, to understand the global structure 
of the halo, sec Section 14.2.1 1 Note that this dataset is not optimal to search 
for substructure in the Solar neighborhood because of the relatively large 
velocity errors, which are typically in the range 30-50 km s~^ ). 

So far, there have only been a few detections of substructure near the Sun. 
This is not really surprising. Because of the short dynamical timescales in 
the inner Galaxy (a typical halo star near the Sun will have made between 50 
and 100 revolutions in a Hubble time), each accreted satellite will generally 
give rise to multiple streams (or wraps of debris) that will cross the Solar 
neighborhood, as shown in Figure [5] In a stellar halo completely built from 
disrupted galaxies, the predicted total number of nearby streams is 

reams 

{300 - 500}(i/10Gyr)3 (9) 

i.e., a few hundred st reams are expected t o cross the solar neighborhood 
after a Hubble time ( Helmi fc White! Il999f l. This quantity is only weakly 



dependent on the number and mass spectrum of the accreted satellites. Fur- 
thermore, in the context of the ylCDM model, it is worthwhile stressing 
that these streams are expecte d to originate in only a handful of objects 
(iHelmi. White fc SDringe"ill2n03f ). 

To detect these streams, samples with approximately 5000 halo stars are 
needed (that is, A'sampic ~ 10 x iVstreams)- Given the size of the samples 
currently available, we could clearly be in the limit in which essentially each 
stream is popula ted by just one sta r. 

Nevertheless, I Helmi et al.l (|l999^ successfully detect ed two strean i s con i- 
ing from a common satellite accreted 6 — 9 Gyr ago ( Kepley et al.ll2007| l. 



A simulation to reproduce its properties is shown in Figure \5[ The reason 
that these authors were able to detect the streams is three-fold. Firstly, they 
studied the distribution of stars in an "integrals of motion" space (defined 
by two components of the angular mo mentum), where stars from a common 
progenitor are very strongly clustered ( Helmi fc de Zeeuwll2000[ ). Thus they 



could map the two streams, one with 4 stars and one with 8 stars, to one 
overdensity with 12 stars. The second reason for the successful detection is 
that the stream stars are going through their pericenters, and at the orbital 
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turning points one expects an enhancement in the density of a stream. (This 
phenom enon is in es sence the same as the shells observed aro und elliptical 
galaxi es jQuinn|[l98^ . This explains why the streams detected bv lHelmi et alj 
(|1999D were not really that cold. The third reason is that the values of the 
integrals of motion of the accreted stars are quite different from those typical 
of solar neighborhood halo stars -lying in the tails of the distribution, exactly 
wh ere one expects the m ost recently accreted objects to be. 

ISeabroke et aL I (|2008 f) have recently put constraints on the presence of 
nearby massive halo streams using the radial velocities of stars towards the 
Galactic poles. They find no vertical flows with densities greater than 1.4 
stars/kpc'^. This rules out the possibility that streams from the Sagittar- 
ius dwarf are crossing the Solar neighborhood, which in tur n implies that 
models in which the dark matter halo is prolate are favored ( Helmil [2004bl : 
iLaw. Johnston fc Maiewskill2005t l. 

Se abroke's results are in excellent agreement with the arguments of lGouldl 
who constrained the degree of granularity using the revised New 
Luyten Two- Tenths (rNLTT) catalog of proper motions of 4500 nearby halo 
stars. He finds, using statistical arguments, that no streams crossing the Sun 
contain more than 5% of the nearby halo stars. This is also in unison with 
the results of Helmi et al. (1999). 

As mentioned above, further detections of substructure will require signif- 
icantly larger samples with full and accurate 3D kinematics. The first steps 
towards this will come with RAVE and SEGUE, but the cle ar breakthrough 
will b e made by the European Space Agency's mission Gaia ( Perrvman et al.l 
Gaia is scheduled for launch in 2011 and will provide high accuracy 
astrometry (at the 10 microarcsecond level for V < 12 stars, and 200 mi- 
croarcseconds for V ^ 20) and multi-color photometry for over 1 billion stars 
throughout the Milky Way. Radial velocities at 10 km accuracy will be 
acquired for stars down to 1^ ^ 16 — 17. Therefore, the Gaia dataset will 
allow us to quantify the amount of substructure in the halo of our Galaxy, 
and consequently its evolutionary path. 



3.3 Implications 

Much of the work in the field of "substructure" has so far been rather qual- 
itative in nature. The discoveries of the substructures are clearly the first 
step, but what do they actually tell us about how our Galaxy formed? This 
re quires a mo re quantitative approach. 

et al.l (j20I)8) have used SDSS data to compute the rms fluctua- 
tions in the number of objects compared to a smooth halo model. They 
find typical values of ^ 4 0%. They interpret this using the simulations of 
iBullock fc .TohnstonI (I2005D . where this statistic can vary from 20% to 80% 



depending on the detailed halo formation history. This could imply that the 
(outer) Galactic halo is fairly average in its history and could have been 
built completely by accreted satellites. However, a careful comparison of the 
overall structure of the fluctuations shows that there may still be room for 
a smooth halo component, as illustrated in Figure [H For example, there are 
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Model 9 Qist 



Fig. 8 Residuals (SDSS or simulations minus the smooth model) smoothed using a 
a = 42' Gaussian from the best oblate model fits for the SD SS data (top left panel) 
and for the 11 simulations from Bullock fc Johnstoiil (120051 ). Only the 20 < r < 20.5 
slice, corresponding to heliocentric distances ~ 14kpc are shown. [From lBell et al.l 
(f200a ). Courtesy of Eric Bell]. 



significant under-densities in the lBuUock fc JohnstonI ( 2005( ) maps, which are 
not apparent to the same extent in the SDSS data. 

IStarkenburg et al.l (|2008D use a different approach and data from the 
Spaghetti Survey to quantify the amount of substructure. Their sample con- 
tains approximately 100 red giant stars with radial velocity and distance 
information. They develop a statistic called the 4-distance, which is essen- 
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tially a correlation function in a 4D space: 



4dy = ^ a^(cjy )2 + ad{di - djY + a„(Urad,i - fradj)^ (10) 

where denotes the angular separation on the sky of star i with respect 
to star j, di its distance and Urad,i its radial (line-of-sight) velocity. The 
coefficients denote weights to normalize the various observables and to 
account for the observational errors. Substructure would manifest itself in 
an excess of nearby pair s in this 4 dimensional space. Using this method, 
[Starkenb urg et all ' (|2008t ) find that - 20% of the stars are in substructures, 



the mos t prominent of which corre sponds to the Sgr stream northern over- 
density (jPohm-Palmer et al.ll200ll ). When comparing to the simulations of 
halos b uilt up completely from a ccreted satellites of 10'' Mq bv lHarding et al.1 
()200l[ ). [starkenburg et al.l ( 2008f ) find that these typically show a significantly 



larger amount of substructure. Only in about 10% of the cases, the character- 
istics of such halos are similar to what is measured by the Spaghetti survey. 
They interpret this as meaning that if the whole halo was built from accreted 
satellites, these must have left behind relatively broad streams, implying that 
either they were relatively large or have orbited the Galaxy long enough to 
be significantly broadened. The other possibl e interpretation is that th ere is 
an underlying smooth component in the halo. lStarkenburg et al. I (I2008D were 



clearly limited by the small sample size, which provided significant Pois- 
son noise to their statistical estimations. This situation is likely to change 
in the very near future with SEGUE, which will provide orders of magni- 
tude larger samples of red giant stars with accurate distances and velocities 
(fMorrison et al...2008b ). 



4 Ages and chemical abundances of halo stars 

4.1 Age distribution 

The determination of the age of individual stars is one of the most difficult 
measurements in astrophysics. Some possible techniques include isochrone fit- 
ting of individual stars (their absolute magnitude, effective temperature and 
metallicity needs to be accurately known) or star clusters; asteroseismology 
and nucleochronology. 

Direct age measurements can be inferred by comparing the observed 
abundance ratios of the radioactive elements of e.g. ^^^Th (half-life of 14 
Gyr) and ^38^ (4.5 Gyr) with theoretical predictions of their initial produc- 
tion. The first measure ments of very meta l-poor nearby stars have yielded 
ages of 15 2 ± 3 .7 Gyr (ISneden et al.lll99fih. 14 ± 3 Gyr (jGavrel et al.ll200ll : 



iHill et al.l[200l . and 13.2 Gvr (|Frebel et al.i r20071 



The white dwarf luminosity function, and especially its cut-off, also pro- 
vides str ong constraints o n the age of the population to which they belong 
(see, e.g. IChabovCTlll998[ ). This method, based on the white dwarf cooling 
sequence, has been used i n the past to derive the age of the oldest stars in 
the Galactic disk (see e.g. iHansen fc Liebertl[2003h . However, the scarcity of 
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known halo white dwarfs and the lack of good kinematical data to distin- 
guish these from the far mor e common disk star s have so far, prevented such 
a measurement for the halo ( Torres et al.111998! ). 

The horizontal branch morphology of a stellar population at a given abun- 
dance tends to be bluer for older populations, and redder for younger ones 
(|LavdenllT998bh . In the halo, there is some indication that the horizontal 
branch morphology of field stars becomes redder with Galactocentric dis- 
tance. If this is interpreted as due to age, it means that the ou ter halo would 
be younger by a few Gyr, compared to the inner halo (see e.g. lSearle fc Zimil 
119781 : ISunt zeff et al]ll99lb . It should be noted, however, that the existence of 
an age spread or a gradient in the population of halo globular c lusters is still 
under debate (e.g. iSaraiedini et al...l997. : .Stetson ct al...l996;.De Angeli et al] 
1200,51 ). 



4.2 Chemical abundances 
4.. 2.1 Metallicity distribution 

The local halo metallicity distribution peaks at a value [Fe/H] — 1.6 dex, 
and extends well below [Fe/H] ^ —3 dex (e.g. iRvan &: Norr"y |l996f). 

Recently ICaroUo et al. based on SDSS data revived the idea that 

the halo may be described by two broadly overlapping components. The inner 
halo peaks at [Fe/H] ~ —1.6 dex, is flattened and slightly prograde, while 
the outer halo would peak around [Fe/H] ^ —2.2 dex, would be rounder and 
in net retrograde rotation. Although a dichotomy in the Galactic halo had 
been suggested before (see also Section [^TT|) . this has now been confirmed 
with a much larger sample of stars covering a large area on the sky and is 
now clearly reflected also in the metallicity distribution as shown in Figure [9l 

Because of the importance of metal-poor stars in constraining the for- 
mation history of the Galaxy, as well as the physics of the high-redshift 
Universe (e.g. reionization, initial mass function, etc) significant efforts have 
been made to better characterize the metallicity distribution at the very 
metal-poor end. More recently, the HK and the Hamburg/ESO (HES) sur- 
veys have yie lded a dramatic increase in the number of nearby metal-poor 
stars known ( Beers fc Christliebl l2005f ) . In these surveys, metal-poor can- 



didates are selected from objective-prism spectra for which the Ca H and 
K lines are weaker than expected at a given {B — V) color. These stars 
are then followed up with medium-resolution spectroscopy (see, for example, 
lEebeTeTaDHQOi). 

These surveys have demonstrated two important points. First, three ex- 
tremely metal-poor stars with [Fe/H] < —4.75 dex have been discovered, all 
of which are carbon rich. Secondly, no metal-free stars (Population III) have 
been found thus far. 

These results could in principl e, put strong constrai nts on the initial 
mass function (IMF). For example. ISalvadori et al. I (|2007D find, using semi- 



analytic models of galaxy formation, that the lack of metal-free stars con- 
strains the critical me tallicity for low-mass star formation to be Zci- > 0, (see 
ISchneider et al.ir2002[ l. or the masses of the first stars to be Mpopni > 0.9 M©- 
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iBrook et al.l (jloo3) reach a similar conclusion using cosmological liydrody- 



namical simulations. 

From a theoretical point of view, the initial mass function of Population 
III stars is expected to b e top-heavy, with m asses in the range 100—1000 Mq, 
e.g. lBromm et al ] (ll999ll : lAber et al Such stars are predicted to leave 



very characteristic imprints in the abundan ce patterns of the next ge neration 
of stars (a very strong odd-even effect, e.g. iHeeer fc WooslevI (l2002D 'l. How- 



ever, such pa tterns are not consistent with those observed for example by 
ICavrel et all (|2004) for nearby very metal-poor stars. Furthermore, the ob- 
served trends can be explained by normal supe rnovae with masses < 100 M0 
(jUmeda fc Nomotdl2005l: ISilk fc Langeill2006D . 



The solution to this conundrum is still to be found. Given the small num- 
ber statistics, the IMF appea rs to be consistent wit h Salpeter as well as with 
a top-heavy functional form ( Salvador! et al.ll2007t ). On the other hand, the 



high incidence of carbon-enhanced, unevolved stars among extremely metal- 
poor stars has been used to p ropose that the IMF in the early Galaxy was 
shifted toward higher masses ( Lucatello et al.ir2005t ). Based on the idea that 



the two most metal-poor stars known could be both in binary systems (and 
that their high [C/Fe] could be due to the winds from an AGB companion 
star), Tumlinson (2007) discusses that an IMF with a characteristic mass 
^ 3 Mq would reproduce the properties and distribution of the halo stars at 
the metal-poor end. 

In the coming years, our knowledge of the halo metallicity distribution 
function is likely to improve significantly with surveys such as SkvmappeiF^. 
SEGUE and LAMOStEB, which will extend b y large factors t he currently 
known handful of extremely metal-poor stars ( Christliebl l2006| ) . These new 



datasets will give us to access to the very high-redshift Universe through the 
fossil record left by the first generations of stars ever formed. 

^.2.2 Elemental abundances 

Much more information about the chemical history of the Galactic halo may 
be obtained from detailed studies of the elemental abundances patterns. This 
is because the various chemical elements are synthesized by di fferent pro- 
cesses in stars of different masses and on different timescales (jMcWilhaml 



For example, a-elements (O, Mg, Si, Ca, S, and Ti) are produced 
during the explosion of a massive star as a Supernova Type II, which occurs 
only a few million years after its formation. Iron-peak elements, on the other 
hand, are produced both by Type la and Type II supernovae. Type la su- 
pernovae are the result of a thermonuclear explosion induced by the transfer 
of mass onto the surface of a white dwarf by a binary companion star. This 
implies that these explosions take place typic ally on a longer timescale, of the 
order of 0.1 to a few Gyr (e.g. iMatteucci fc Recchill200ll ). Heavier elements 
beyond the iron-peak are created by neutron capture, through slow (s) and 
rapid (r) processes. The s-process can take place when the neutron flux is 
relatively low, i.e. the timescale between neutron captures is large compared 

http:/ /www. mso.anu.edu. a u/skymapper/l 
(http://www.lamost.org, 
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Fig. 9 Metallicity distribution function based on medium-resolution spectroscopy 
for 2 X 10'' nearby stars obtained with SDSS. The histograms show that the inner 
halo (V ~ —200 km s~' ) exhibits a peak metallicity [Fe/H] — —1.6, and that the 
outer halo (right panels) peak at a value [Fe/H] = —2.2, and have a net retrograde 
rotation. [From CaroUo ct al. ( 2007,) . Courtesy of Daniela CaroUo. Reproduced with 
permission of Nature.] 



to that of the beta decay. These conditions are found in the envelopes of 
Asymptotic Giant Branch (AGB) stars, and are most efficient in stars with 
mas ses 3 — 5M(T), iniplying a timescale of 100 Myr after the stars were born 
(see iGorielv fc Siesi 120051 and references therein). On the other hand, the 
r-process occurs when the neutron flux is sufficiently high to allow for rapid 
neutron captures. This is believed to occur in environments like those pro- 
duced by supernovae Type II. It should be stressed, however, that the exact 
sites and conditions under which the various neutron-capture elements are 
produced, particularly at very low metallic ities, are still somewhat uncertain 
(|Travagho erani2004 lArnould et al.ll2007t) . 

Figure [10] shows how the chemical enrichment depicted by the [a/Fe] 
versus [Fe/H] may proceed in a given object depending on its initial mass 
function (which sets the initial level of [a/Fe]) and its star formation rate 
(which sets at what metallicity the "knee" is observed). 
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Fig. 10 Sketch of the predicted [0/Fe] versus [Fe/H] relation in different systems 
as a consequence of their different enrichment timescales. The a enhancement seen 
in the Galactic bulge can be explained by an intense star formation burst in which 
mainly supernovae Type II polluted the interstellar medium. On the other hand, the 
Magellanic Clouds had a lower star formation rate, and so most of their stars have a 
relatively enhanced iron-content. Thin-disk stars near the Sun are an intermediate 
case. [Based on .Matteucci, (,2001. ').] 



Interestingly, the majority of the local halo stars are a-enhanced ([Wheeler et alJ 
119891: iNissen et al.lll994 ICarretta et al.ll2000t) . where typically [a/Fe] ~ 0.3 
dex. This implies that the formation timescale for these stars was very short, 
less than 1 Gyr. This argument is based on the fact that the a-enhancement 
may also be seen as an iron deficiency. Such a configuration may be obtained 
when the interstellar medium (ISM) has been enriched by s upernovae Type 
II, but not significantly by Type la (jGilmore fc Wvselll998| ). 

More detai l ed ins ights have been provided bv ICayrel et al.l ( 20041 ) and 
[Francois et al.l ( 20071 ). They have obtained high- resolution sp ectra (R ^ 
45000) for 30 very metal-poor giants (-4.1 < [Fc/H] < -2.7 dex). [Cavrel et al.l 
(|2004[ ) find that the scatter in several [a/Fe] and [Fe-peak/Fe] ratios is con- 
sistent wi th being only due to observational errors, as shown in Figure 1111 
see also ICohen et al.ri2004l) . This would suggest that these stars formed 
in an extremely well-mixed medium, i.e. the cosmic scatter expected from 
the pollution by supernovae of different masses is not observed. In contrast, 
the neutron-capture elements behave ve ry differen tly, showing a large spread 
at low metallicity ( Barklcm et al. 20051 iFranpois et al. 2 007, , and references 
therein). It is currently unclear how to reconcile these results. In a scenario 
in which the halo was formed via the accretion of several independent ob- 
jects, one would naively expect that each of these objects followed their own 
chemical evolution path, hence a superposition of their debris should give rise 
to a large scatter in the various chemical abundance patterns. Perhaps the 
observed large scatter in [neutron-c apture/Fe] would be consistent with such 
a scenario (see ICescutti et al ] |2006l ). but why is there not a spread in other 
elements? Uncertainties in how nucleosynthesis proceeds at very low metal- 
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Fig. 11 Elem ents ratio s as function of [Fe/H] for the sample of very metal-poor 
stars from Ca vrel et al.l (2004). Ca is an a element, while Cr is from the iron- 
peak, yet in both cases the scatter is extremely small. In contrast, neutron-capture 
eleme nt ratios such as [Y /Fe] a nd [Ba/Fe] show ve ry large scatter in the same 
stars ([Francois et all 120071 ). [From [Cavrel et all (|2004l ). Courtesy of Rogier Cayrel. 
Reproduced with permission of A&A.] 



licities, as well as in modeling the chemical evolution in the proto-galaxy in 
a cosmological context, are likely to be both elements of this puzzle. 

The chemi stry of stars can also be used to trac e the merging history of the 

Galac tic halo ( Freeman fc Bland-Hawthornl2003 ). For example . iNissen fc Schust"ei1 



(|1997[ ) found a correlation for local halo stars such that those with large apoc- 
entric distances have relatively low [a/Fe]. This may be an indication that 
these anomalous halo stars have been accr eted from dwarf galaxies with a dif - 
ferent chemical evolution history (see also lFulbrightl[2003 llvans et al.ll2003t ). 
This is consistent with the fact that stars in dwarf spheroidal galaxies, tend 
to have lower [a/Fe] at a given metallicity compared to the halo. This in turn 
has been used to show that dwarf galaxies could not have been significant 
contributors of the st ellar halo if accreted after they formed a significant 
fraction of their stars (jShetrone et al.ll200ll: iTolstov et al]|2003t IVenn et all 



[2004). 

More recent abundance studies of large samples of stars in several dwarf 
spheroidals have demonstrated that these objects have similar initial levels of 
[g/Fc ] at very low [Fc/H], as shown in Figure [12] (Tolstoy et al. 2006; LetartJ 
120071 : IhIU et al.il2008f) . It is hkely that their relatively low star-formation 
rates have led to slow chemical enrichment, implying low values of [a/Fe] 
compared to the h alo for all but the ol dest (or most metal-poor) stars in the 
dwarf spheroidals ( Tolstov et al.l[2003| l. 



5 The satellites: dwarf galaxies and globular clusters 

Approximately 150 globular clusters and about 20 satellite galaxies are found 
in the Galactic halo. It is interesting to study their distribution as well as 
their internal properties to elucidate the evolutionary path of the Galaxy. In 
particular, it is worthwhile comparing these objects' stellar populations to 
the halo field stars, since in the most naive interpretation of the hierarchical 
paradigm, the stars have their origin in small galaxies, presumably akin to a 
certain (to be defined) extent the dwarf satellites we see around us today. 
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Fig. 12 Q-abundance (average of Ca, Mg and Ti) as functi on of [Fe /H] for 92 
member stars in a central field of Sculptor (crosses, from Hill ct al."2008') compared 
to those in our Galaxy (dots; from the compilation of Venn ct al. 2004). The open 
squares are previous measurements of individual stars tShet rone et al. 2003). The 
most metal-poor stars in the dwarf spheroidals are as a enhanced as field stars in 
the Galactic halo. It is only at higher [Fe/H] that the [a/Fe] is significantly lower in 
the dwarf spheroidals (jToIstov et al.j|2006i ). [Courtesy of Eline Tolstoy and Vanessa 
Hill.] 



The number of known dwarf galaxies in the Local Group has doubled 
in the past few years, particularly at the faint luminosity end, thanks to 
the discovery of a large number of very low surface brightness satellites (e.g. 
iBelokurov et al.|[2007at iMartin et al.ll2007l ). Despite this recent large increase 
in the abundance of Galactic satellites, this is still well-below the number 
of bound su bstructures orbiting g alaxy-size halos in pure cold dark matter 
simulations ( Koposov et al.|[2007). This discrepancy has been ter med the 
"missing satellites problem" ( Klvpin et al.lll999t iMoore et al.lll999t) . Several 
solutions have been put forward including, for example, a modification of the 
nature of dark matter (Spcrgcl fc Stcinhardt 2000 ). On the other hand, the 
modeling of the baryonic physics at the low-mass end of the galaxy scale is far 
from complete, and this is where environmental effects and processes such as 
reionization, feedback from massive stars, etc. are likely to have a significant 
impa .ct on the probability of a substructure to hold onto its gas and form s tars 
(e.g. iMac Low fc Ferraralfl999l : iBullock et ani2000l : iKravtsov etallbOO^ ). 

The newly discovered satellites have luminosities in between (although 
sometimes comparable to) globular clusters and the classic dwarf spheroidals, 
but they are significantly more ex tended than the form er. They have ex- 
tremely large mass-to-light ratios (|Simon fc Gehal [20071 ). computed under 
the assumption of virial equilibrium (which may very well turn out to be a 
very poor approximation since these objects could well be tidal debris). 

There are two very interesting properties of the spatial distribution of the 
satellites. Firstly, their distribution is rather centrally concentrated, in the 
sense that th e median distance is less th an 30% - 40% of the virial radius (as 
computed bv lBattaglia et al.ll2005Ll2006l ). For comparison, satellite substruc- 
tures around halos in dark matter cosmological simulations tend to be much 



28 



Amina Helmi 



more evenly distributed (jGao et alJ 120041 : iKang et all l2005l ) . This apparent 
discrepancy can be solved by assuming, for example, that the luminous satel- 
lites are embedded in those substructures that were most massive at the time 
they crossed the virial radius of the Galaxy. Also the radial distribution of 
the halo globular clusters follows very closely that of the field stars. This 
can be explained if these obje cts formed in a highly bia sed set of gas-rich 
proto-galaxies at high-redshift ( Kravtsov fc Gnedinll2005l ). 



The second in teresting property is that the satellites are quite anisotropi- 

cally distributed (|Hartwickl2000HKroupa. Theis fc Boilvl2005HM"etz. Kroupa fc Jerjenl 



l2007f ). This can be reconciled with cold dark matter models if these galaxies 
have fallen in from a fe w directio ns, perh aps in gro ups (Li & Hclmi 20081), or 
along filaments (Libeskind et al.L 2005; Zentner et al. 2005; Kan g et al. 2003 ). 
It is unclear whether this peculiarity will remain once the sky coverage is 
more complete, and once a thorough survey of the faintest objects is avail- 
able. Further more, establ ishing the reality of the dynamical families proposed 
by Lvnden-Bell &: Lvnden-Belli (|l995l ) , is only possible with knowledge of the 
satellite's orbits. For this, precise tangential velocities (i.e. proper motions) 
are required. Unfortunately only a few satellite galaxies have proper mo- 
tion measurements, mostly obtained using the Hubble Space Telescope (see 
e.g. iPiatek et al . 2007i , and references therein) and these are still affected by 
large uncertainties. However, in the coming decade, the situation is going to 
improve significantly thanks to Gaia. 

In contrast, the dynamics of the (in ner) halo globular c luster popula- 
tion is much better known. For example, iDinescu et al has compiled 



a sample with > 40 clusters with well-measured tangential velocities. The 
measurements are relatively modern (i.e. made with high-precision measur- 
ing machines, such as used in the Southern Proper Motion survey), and are 
tied to an absolute reference frame using, e.g. galaxies and QSOs, or stars 
with known absolute proper motion from, e.g. the Lick No rthern Proper Mo- 
tion Program or Hipparcos (for more details, see e.g. IDinescu et al.lll999l . and 
references therein). 

The reliable measurements of the dynamics of globular clusters in com- 
bination with dedicated programs to characterize their properties, such as 
chemical abundance patterns, ages, horizontal-branch morphology, and struc- 
tural parameters has prom pted the search for correlations in the hope to 
trace Galactic history (e.g. Ashman fc Zepf '1998: 'M ackev fc Gilmorell2004L 
a nd references therein) . Thi s is much in the same spirit of the original work 
of ISearle fc Zl^ (|1978D andlZinn| ([l99l). 



A recent highlight is the dis covery of multip le stellar populations in the 
most massive globular clust ers (iLee et al. 19991). O ne proposed explanation 



(but not the only one, c.f. iBaumgardt et alT 20081 and references therein) 



is that these might be remaining cores of disrupted nucleated dwarf galax- 
ies, as is the case for M54 which is located in the Sagittarius dwarf itself 
( Layden fc Saraiedinill200dh . These pecuhar globular clusters also show an 



extended horizo ntal branch (EHB) , and the presence of double or broadened 
main sequences ( Piotto et al.ll2007t ). In fact, about 25% of the globular clus- 



ters in our Galaxy exhibit unusual EHBs, which has prompted the idea that 
they might have a unique origin. Furthermore, these globular clusters are 
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Fig. 13 Star formation histories of Milky Way dwarf spheroidal satellites. Each 
population box gives a schematic repr esentation of s tar formation rate (SFR) as a 
function of age and metaiiicity. [From [Grebell (|2000l ). Courtesy of Eva Grebel, and 
Fabio Favata (editor).] 



on average more massive and share the ra ndom kinematic s characteristic of 
normal globular clusters in the outer halo (|Lee et al.ll2007f ). 

These results, as well as the evidence of dyna mical relation between the 
Sagit tarius dwarf and six Galactic globular clusters ( Bellazzini. Ferraro fc Ibatal 



|2003[ ). suggest an acc retion origin for the outer halo clusters. This is also sup- 
ported by the work of lMackev &: Gilmord ( 2004t ). who conclude that all young 



halo clusters (with [Fe/H] < —0.8 dex and red HB morphology, more com- 
mon in the outer halo) and ~ 15% of the old halo clusters ([Fe/H] < —0.8 
dex and blue HB morphology) are of external origin. 

Dwarf galaxies show, of course, much more complex stellar populations. 
Furthermore, these galaxies depict a larg e variety in star for mation and 



i^urthermore, triese galaxies depict a. larg e variety m sta.r for mation and 
chemical evolution histories ( Grebell I2OOOI : iDolphin et al.|[2b05() . as shown 



in Figure [THl Some dwarf galaxies appear to have stopped forming stars very 
early on, although it is clearly difficult to disentangle whether this happened 
at z ^ 10, i.e. before reio nization, or at z ^ 3 — 5, i.e. after reionization 
(|Gnedin k Kravtsovl [20061 ) . On the other hand, some dwarf gala xies have 
mana ged to form st ars until very rec ently, such as Fornax (jStetson et al.l 
|1998[ ) or Sagittarius (jSiegel et al.ll2007[ ). 



There is, however, a remarkable uniformity in that all dwarf spheroidals 
studied so far show an absence of sta rs with metallicities below a value of 
[Fc/H] 3 dex (|Helmi et al.l I2006D . This could suggest a uniform pre- 



enrichment across the Local Group at high redshift, before the first genera- 
tions of stars in these objects formed. 
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We are also interested in understanding the link between the f ield stars in 
the h alo and those still bound to the satellites of the Milky Way ( Sales et all 
I2007D . For example, the contribution of present-day dwarf galaxies to the 
halo near the Sun can be ruled out on the basis of the colors and chemica l 
properties of the stellar populations ( Unavane et al.lfl996l : IVenn et al.ll2004^ ■ 
These dwarfs have had a Hubble time to evolve in relative isolation, and it is 
reasonable that their stellar populations would not resemble those of objects 
disrupted a long time ago. However, one may wonder what the relation is 
between the progenitors of the dwarf galaxies and the halo building blocks. 
Under the assumption that the very metal poor stars represent the first 
generations of stars formed, the apparent lack of such objects in the dwarf 
spheroidals compared to the halo, implies that their progenitors were funda- 
mentally different. One possible scenario is that the halo building blocks were 
associated to high-density peaks that collapsed at higher redshifts, while the 
dwarf galaxies presu mably descend from more average densit y fluctuations 
in the early Universe (|Helmi et al.l 120061 : ISalvadori et al.ll2008[) . 

Nevertheless we know at least of one example of a significant contributor 
to the outer halo: the Sagittarius dwarf galaxy. Although it is currently un- 
clear how massive this object was initially, current models postulate an initial 
stella r mass of ~ 2 — 5 x 10* M© (Helmi 200_4a.; Law. John ston & Maiewskl 
which is slightly smaller but comparable to the total mass of the stellar 
halo obtained by extrapolating the local density and assuming a power-law 
distribution as described in Section [2l Therefore, given that we are just be- 
ginning to map the outer stellar halo in sufficient detail, it is very likely 
that much of this region of the Galaxy is the result of the accretion of dwarf 
galaxies like those we observe today. 



6 Formation scenarios 

There has been a tendency in the literature to classify formation scenar- 
ios according to two traditional models introduced in the 1960s and 1970s 

respectively, to describe the evolutio n of the Galaxy. 

The first model was proposed bv'Eg gen. Lvnden-Bell fc Sandaed ( 1962I) . 



where "the oldest stars were formed out of gas falling toward the galactic 
center in the radial direction and collapsing from the halo onto the plane. 
The collapse was very rapid and only a few times 10* years were required for 
the gas to attain circular orbits in equilibrium" . The influence of this work 
has been phenomenal, not just on models of the stellar halo, but also in the 
field of galaxy evolution in ge neral. The second mo del, which is often used as 
a counter scenario, is that bv lSearle fc Zi^ (|1978D . On the basis of the lack 



of an abundance gradient for the outer halo globular clusters, these authors 
proposed that "these clusters formed in a number of small protogalaxies that 
subsequently merged to form the present Galactic halo" . For many authors, 
this latter model is more reminiscent of the currently popular hierarchical 
paradigm. However, as we shall see below, ingredients of both models are 
likely to be present in the assembly of a Galactic system in a hierarchical 
cosmology. 
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It is therefore perhaps more instructive to distinguish models according 
to the degree of diss ipation involved in the formation of a given component. 
In this context, the lEggen. Lvnden-Bell fc Sandagel ( 19621 ) model would fit 
in a category where the large majority of the nearby halo stars were formed 
during a dissipative process (name l y the collapse of a gaseous cloud). On the 
other hand, in the ISearle fc ZinnI (|l978h model, the assembly of the outer 
halo would have proceeded without any dissipation. 



6.1 Global properties 

Several authors have attempted to model the formation and characteristics of 
the Galactic stellar halo. Different techniques have been applied, such as nu- 
merical simulations with and without gas physics and star formation, which 
are sometimes (fully) embedded in a cosmological context. Other attempts 
include the use of phenomenological descriptions of the evolution of baryons 
inside halos, usually in combination with an N-body simulation that pro- 
vides the dynamical history of the system in the context of the hierarchical 
paradigm. 

The global properties of the Galactic stellar halo, such as its shape, cor- 
relations with kinematics and metallicity distribution, have been reproduced 
in many of these works. 

For example, ISamland fc GerhardI (|2003D have modeled the formation of 



a disk galaxy in a slowly growing dark matter halo. In essence, this model is 
quite similar to th cEggen. L vnden-Bell fc Sandage (1962.) scenario, but takes 
into consideration the fact that galaxies have grown in mass from the Big 
Bang until the present time. They follow the cooling of gas, star formation 
and feedback processes in great detail, and are able to obtain a realistic disk 
galaxy at the present time. In this simulation, the stellar halo defined by 
the stars with [Fe/H]< —1.9 dex is not rotating, while the more metal-rich 
stars with -1.9 < [Fc/H] < -0.85 are rotating at ~ 70 km s^^ . The first 
component is slightly flattened in the same sense as the disk, while the second 
component is significantly flatter and largely supported by rotation. In this 
sense, one may argue that it resembles more closely a thick disk than the 
inner halo of the Milky Way. 

More realistic attempts to simulate the formation of a disk galaxy like the 
Milky Way in a cosmological f rame work include t hose o f fSteinmetz fc Miilleil 
(|1995[ ): iBekki fc Chibal (|200lD and iBrook etall (|2003[ ). These authors have 



modeled the collapse of an isolated, rigidly rotating sphere of 10 M0 
onto which small-scale fluctuations according to a cold dark matter spec- 
trum have been imposed. In these works the evolution of the dark matter 
and of the gas, star formation, etc. are followed self-consistently. In this 
sense they are, in principle, able to determine the relative importance of dis- 
sipative processes versus dissipationless mergers in the build-up of the halo. 
However, the numerical resolution in these studies is low (the spatial resolu- 
tion, as measured by the softening, i s typically "^ 2.5 kp c), and the number 
of particles is small. For example, in iBrook et al. 1 dool) the dominant disk 



component and the spheroid at the final time are represented by ~ 4 x 10^ 
star particles. This implies that the resolution is not sufficient to determine 
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the detailed properties of the simulated galaxies robustly (e.g. small scale 
substructure such as dwarf galaxies and streams are not well resolved) . Nev- 
ertheless, global (average) properties are likely to be reliable. For example, 
ISteinmetz fc Miillej (^1995^ find in their simulations that the stellar halo (de- 
fined as the more metal-poor component) is triaxial in shape (c/a = 0.65 and 
h/a = 0.85, with the minor axis roughly aligned with the angular momen- 
tum of the disk), and th at the density profile is well- fit by a de Vaucouleurs 
law. In the simulation of iBekki fc Chibal (|200l D the outer halo is completely 
built up through accretion, while the inner halo results from merging be- 
tween gas rich clumps. T he inner halo is also adiabatically contra cted by 
the formation of the disk (|Binnev fc Mavlll986l : IChiba fc Beei^l200lD. which 
explai n s its more-flattened s hape. In agreement with ISteinmetz fc Miillerl 
(|1995( ). iBekki fc Chibal (|200lh find that the density profile varies from r''^ 
at r < 10 kpc to r~'^ at larger distances. 

This is also in agreement with the simulations of lAbadi et al. I dlool). 
These authors have run fully cosmological hydrodynamical simulations of the 
formation of a disk galaxy. Although their final galaxies are of too early Hub- 
ble type (most of them have bulges that are too dominant compared to the 
disk component), they find some common (and hence presumably general) 
characteristics. The surface brightness profile of the outer halo components 
of these galaxies (r > 20 kpc) can be fit with a Sersic law, with (n) = 6.3 
and (i?cff) = 7.7 kpc, i.e. R~^-^ at ~ 20 kpc and R^^-^ around the virial 
radius. The shape of the outer halos is very mildly triaxial {b/a) — 0.91 and 
(c/a) — 0.84. These outer regions (r > 20 kpc) contain mostly accreted stars 
(with fractions greater than 90%), which originate in the most massive merg- 
ers that the central galaxy has experienced over its lifetime (the contribution 
of small accretion events may be underestimated because of resolution, but 
this is nevertheless unlikely to be the major channel by which galaxies grow). 
The inner halos ar e formed in-situ and are much more flattened c/a ^ 0.62 
(|Abadi et al.ll200(ih . 

iDiemand et aL I (|2005[ ) provide an alterna t ive vi ew from dark matter only 
cosmological simulations. Following iMoord (|200lh . they identify the most 
bound particles in halos that have collapsed at very high redshift (i.e. 2.5tT 
peaks at z = 12) and find that these end up in the inner regions of the 
most massive object at the present time. These particles (which are meant 
to represent stars) are distributed according to a very steep density profile r" 
with index n — 3 in analogy to the stars in the stellar halo of our Galaxy. 

A complementary approach has been developed by iBuUock fc JohnstonI 
( 2005[ ) who follow the formation of the stellar halo purely built via accre- 
tion of satellites. The evolution of the baryonic components of the satellites 
(fraction of cold and hot gas, star formation and chemical enrichment his- 
torie s) are foll owed through a phenomenological semi-analytic approach (as 
in White fc Fre nk 1991; KauflP mann et al.. 1993: Somcrville fc Primack 1999). 
The satellites are pla ced on orbits ( consistent with those found in cosmolog- 
ical simulations, e.g. iBenso 3 (|2005M around a slowly growi ng gravitational 
poten tial that models the Galaxy (resembling the approach of lSamland fc GerhardI 
(|2 003 f)). In this sense, they are limited to times after which the central po- 
tential has settled (i.e after major mergers), and so, essentially they are 
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well-placed to tackle the formation of the outer halo, i.e. r > 15 — 20 kpc. 
Through this modeling they find that the density profile of their stellar halos 
vary from p oc at 20 kpc to r""* beyond 50 kpc. 

Th e higher resolution and more detailed modeling of lBullock fc JohnstonI 
( 2005[ ) also allows them to make predictions about t he distribution of c hem- 
ical abundances of the Milky Way outer halo stars. iFont et all (|200fiD find 
that the mean metallicity of the halo depends on the accretion history, and 
that gradients in metallicity of ~ 0.5 dex are not necessarily uncommon in 
their accreted halos. They predict that the stars in the outer halo should have 
lower [a/Fe] than stars in the inner Galaxy. In their models, the most metal- 
poor stars are located i n the inner 10 kpc, w hich at face value is opposite to 
the trend suggested bv lCaroUo et al. ] (120071). 

The hybrid cosmol ogical approach of De Lucia fc Helrnil (12008 ) cornbines 
elements of the work of lDiemand et al.] ( 2005 ) and Bullock fc JohnstonI ( 20051 ) . 
These authors use a high-resolution cosmological simulation of the formation 
of a Milky Way-sized dark matter halo coupled to a semi-analytic galaxy for- 
mation code, to obtain insights into the global evolution of the Galaxy. By 
identifying the most bound particles of substructures in the simulation with 
the luminous cores of satellite galaxies, they are able to follow the structure 
but also the metallicity and age distribution of the "si mulated" stellar halo. 
As shown in Figure [111 also iDe Lucia fc Helrnil (l2008l ) find that the stellar 
halo is much more centrally concentrated than the dark matter halo. Further- 
more, their inner halo is dominated by higher metallicity "stars" , vis-a-vis 
the outer halo, with the transi tion occurring at ^ 1 — 15 kpc from the center, 
reminiscent of the findings o f ICarollo et al 1 (I2007D . 



6.2 Current status 

The discussion above shows that it is possible to reproduce the density profile 
and shap e of the inner stellar halo in a model in which all stars were formed in 
situ as in lSamland fc G erhard' ('2003'), as well as i n a model in wh i ch all stars 
were accreted as in lDiem and et al. (2005) . Bullock fc JohnstonI ( 20051 ) and 



iDe Lucia fc Helml ()2008l ). Furthermore, in all the models discussed here the 
oldest stars are located in the spheroid (halo-|-bulge), which implies that their 
age distribution and the broad chemical properties (low mean metallicity and 
a-enhancement) are also likely to be reproduced. Nevertheless, much more 
work is required to understand the expected trends as well as predictions 
for the detailed chemical abundance patterns that have no w been observed 
i n suc h exquisite detail for nearby halo stars. The work of IChiappini et all 
(Il997h is an excellent example of how this can be done, but it needs to be 



embedded in a more cosmological framework. 

The discussion in the previous paragraphs also shows that we do not re- 
ally possess a clear view of what are the exact predictions of the hierarchical 
paradigm on the scales of the Milky Way and its satellites. The relatively 
simplistic and coarse fully cosmological simulations are limited in their pre- 
dictive power at the level of detail that is required for the Milky Way galaxy. 
The semi-analytic approach of combining high-resolution dark matter simu- 
lations with phenomenological descriptions of the evolution of the baryonic 
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Fig. 14 Evolution of the dark matter distribution in a hiRh -resolution co s moloK - 
ical simulation of the formation of a Milky Way-like halo bv lStoehr et al.l (|2002l ) . 
By identifying the cores of dark matter substructures with luminous satellites, 
IDe Lucia L Helmil (p08) can follow the formation and evolution of a spheroidal 
component, which would correspond to the stellar halo. For example, by analyzing 
where these objects deposit their debris, they can study the link between the chem- 
istry and the structure of the stellar halo. This is shown in the bottom right panel, 
where the green and orange curves denote the surface density profile of metal-poor 
and metal- rich "stars", respectively. In comparison, the dashed black curve cor- 
responds to the dark matter halo, which shows a significantly more shallow and 
extended distribution. [From IDe Lucia fc Helmil POOS ') . Courtesy of Gabriella de 
Lucia.l 
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Fig. 15 Distribution of angular momentum for a sample of nearby halo stars. The 
color-coding indicates different metallicity ranges. Stars with orbits confined to the 

plane (Lperp = ^ L\ -\- L\ < 350 kpc km s~^ ) tend to have a net prograde motion 
(corresponding to Vtf, ~ 50 km s~^ ), and their distribution is very smooth in this 
space. Stars with highly- inclined orbits (Lperp > 350 kpc km s~^ ) have no net 
rotation and their distribution is quite lumpy, indicative of the presence of merger 
debris. It is als o interesting to note th at this "higher halo" is more metal-poor on 
average. [From [Morrison et al.l l|2008al ). Courtesy of Heather Morrison.] 



components, is a powerful and promisin g; technique, that will be developed 
much further in the near future (see e.g. iDe Lucia k Helmill2008[ ). The first 
attempts by iBullock fc JohnstonI ( 20051 ) at modeling the outer halo have 
yielded extremely interesting results, despite the fact that the growth of the 
Galaxy is followed under rather ideal conditions. 

From the observational point of view, there seems to be clear consensus 
that the outer halo is very lumpy. This would imply that there is little room 
for a component predominantly formed via a dissipative collapse. 

In the case of the inner ste llar halo, the situation is much less clear. Re- 
cently, [Morrison et al.l (|2008al ) analyzed the distribution of angular momen- 
tum for a sample of nearby metal-poor halo stars with accurate velocities. 
These authors found a break in the properties of these stars depending on 
their angular momentum as shown in Figure 1151 The distribution of angu- 
lar momentum is clearly not smooth, in particular for the higher halo (stars 
on highly inclined orbits), which suggests some amount of accretion in the 
build-up of this component. 

Perhaps the most direct way of testing the amount of "dry" (i.e. dis- 
sipationless) accretion or mergers that have contributed to the stellar halo 



36 



Amina Helmi 



near the Sun, is to find the debris streams that are expected to be present. 
The granularity in the velocity distribution of stars should be measurable, 
if present, in large catalogs of halo stars (> 5000) with very accurate kine- 
matics (errors < 10 km s""'^). The only processes by which these streams 
are expected to diffuse away are through collisions, or if the gravitational 
potential is (or ha s been) extremely chaotic . Galaxies are essentially col- 
lisionless systems ( Binnev fc Tremaind Il987i) . and the effect of dark mat- 



ter substru cture on the coherenc e of tidal streams has been shown to be 
very small ( Johnston et al.l l2002| ) . Therefore one may only need to worry 



about the chaotic nature of the potential. First studies performed on fully 
cosmological dark matter simulations show evidence that the number of 
stellar streams expected near the Sun is the same as in steady-state mod- 
els of the accretion of satellite galaxies onto a fixed gravitational potential 
(|Helmi. White fc SpringefeoO,'^ . The Galactic bar may also be a source of or- 



bital diffusion. This triaxial pot ential is known to induce resonance s amongst 
stars on low inclination orbits ( Dehnenll200lt [Pamaev et al.ll2005[ ). and may 



well cause streams to become more diffuse in time, and hence harder to re- 
cover. 

Another promising approach is that of chemical tagging, i.e. stars born 
within the same m olecular cloud should have very similar chemical abun- 
dances as shown bv iDe Silva et al.l (j2007f ). We have yet to understand how 
much scatter there can be in a galaxy that was accreted very early on, and in 
particular, if similar degrees of coherence in the abundances are expected for 
such (much larger) systems. For exam ple, dwarf gal axies today show metal- 
licity spreads of the order of 0.5 dex ( Mated fl998l ). Although one expects 



stars in a stream to come from a small region of phase-space, it is unclear 
how different their chemical abundances may be. For example, a satellite 
that has been disrupted only relatively recently will give rise to streams that 
are less chemically homogeneous than one whose streams are dynamically 
older. In the first case, the metallicity spread will likely reflect that of the 
whole galaxy, while in the second case it may reflect only that of a small por- 
tion of the system. In the limit of an infinite amount of time, the stars in a 
given location on a stream should have all been born within the same molec- 
ular cloud, and hence have the same exact abundances, provided diffusion 
processes are unimportant. 

An accurate measurement of the ages of stars in the halo would yield 
unique insights. This would give a very robust test of the degree of dissipation 
involved in the formation of the stellar halo. For example, if the stars were 
formed predominantly during gas-rich mergers, the star formation history 
of the halo would have been largely bursty. In this case, each burst would 
be associated to each one of those gas-rich merger events. If on the other 
hand, the halo formed only via the accretion of a few objects containing 
some stars, each of these objects will have had a different star formation 
history, and hence a larger spread in ages might be expected. Note that in 
the limit in which gas-rich mergers have taken place successively on a very 
short ti mescale, such a model will be e s sentia lly indistinguishable from the 
original [Eggen. Lvnden-Bell fc Sandagd ( 1962f ) proposal for the formation of 



the inner halo, even if the ylCDM model is correct. 
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7 Halos in other galaixies 



Stellar halos are generally such low surface brightness components of galaxies 
that detecting th ese beyond the Milky Way is a major challenge. For example, 
Im orriso estimate the surface brightness of the Galactic halo at the 

solar radius to be Ey ^ 27.7 mag/arcsec^. 

Direct testimony of the difficulty of such an enterprise is given by M31. 
Its spheroid was long believed to be a simple extension of its bulge, since 
"metal-rich" stars were f ound even out to a radius o f 30 kpc followi ng a de 
Vauc ouleurs profile (e.g. iMould fc Kristianl 1986: Br own et al.ll2003l ). How- 
ever, iKalirai et al.l ( 2006( ) and lChapman et al.l (j2006( ) recently discovered a 
metal-poor halo out to 160 kpc that is similar to that of the Galaxy, both in 
terms of its metallicity, radial profile, as well as the amount of substructure. 
Our external p oint of view h as allowed the study of this halo in a much more 
global fashion ( Ibata et al.li2007 ), and has provided spectacular images as 
shown in Figure [161 Just as for the Milky Way, in the halo of M31 we find a 
giant stream as well as a wealth of fainter substructures at surface brightness 
levels of ~ 30 mag/arcsec^. 

M33 seems to have a halo much like that of the Milky Way and M31, 
with a similar metaUicity [Fe/H] — 1.6 dex, despite its smaller total mass 
( McConnachie et al.ll2006l ). There is also some evidence that the Large Mag- 
ellanic Cloud may have a stellar halo ( Minniti et al.l 120031 ). However, the 
origin of this halo is less clear because the LMC appears to have been tidally 
pe rturbed by the Milky Way (its disk is asynimetri c and fairly thick as shown 
by Ivan der Marell l200ll : Ivan der Marel et al.l 20021). S uch strong tides could 
give rise to a halo-like component (|Maver et al. I booih . 

Beyond the Local Group, there have been a few studies of halos using 
surface brightness photometry. This work is very challenging because it re- 
quires characterizing a component several magnitudes below the sky level. 
A few halos have bee n detected in this wa y (e.g. NGC5907, NGC3115, see 
iMorrison et al .l ('l994');'Sackett et al.' (1994)) as well as various tidal features 
(IZheng et aill9 99: Pohlen et al. 2004). A n interesting detection was obtained 
by IZibetti et al.1 (j2004f ) who, stacking images of 1047 edge-on disk galaxies, 
were able to measure a diffuse moderately fiattened stellar halo, well de- 
scribed by a power law density profile. 

An alternative, possibly more powerful approach is to study the resolved 
stellar populations in the outskirts of (disk) galaxies, as frequently done in 
the Local Group. Using the Advanced Camera fo r Surveys (ACS) on HST 
iMouhcine et al.l (|2005b') and lde Jong et al.1 (j2007l ) have been able to resolve 
stars a few magnitudes below the tip of the red giant branch for a sample 
of nearb y disk galaxies (l o cated at distances between 4 and 14 Mpc). For 
example, iMouhci ne et al. I (l2005br ) find that the stellar halos are generally 
old, and typically more metal-poor than the main galaxy. These authors also 
find large variety in their properties. An example is NGC891, often referred 
to as the twin of the Milky Way, whose halo contains a much larger spread in 
stellar populations tha n is apparent for the Milky Way at a similar di stance 
(|Mouhcine et al.ll2007f l. On the other hand, iMouhcine et all (|2005ah have 
suggested that there is a correlation between the stellar halo metallicity and 
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Fig. 16 View of the halo of M31. Red, green and blue show, respectively, stars 
with -0.7 < [Fe/H] < 0, -1.7 < [Fe/H] < -0.7 and -3.0 < [Fe/H] < -1.7. The 
differences in stellar populations between the Giant Stream and the several minor 
axis streams can be seen as striking diffe rences in co l or. At the center there is a 
scaled image of the central of M31. [From lTbata et ahl l|2007l ). Courtesy of Rodrigo 
Ibata. Reproduced with permission of the AAS.] 



the parent galaxy luminosity. Our Galaxy falls off this metallicity-luminosity 
relation by more than 1 de x, i.e. the Galactic st ellar halo would be unusual, 
or too metal-poor (see also lHammer et al]|2007l for a similar statement). 

Some caution is necessary when interpreting these observations, espe- 
cially for the following two reasons. First, in the absence of spectroscopic 
measurements, the chemical properties of the stars cannot be determined re- 
liably because of the well-known age- metallicity degenera cy (a clear example 
of which is the dwarf galaxy Carina. [Koch et all ( 20061 )). Secondly, stellar 



halos can be very lumpy, particularly at very low surface brightness levels, 
as shown in Figures Bl and [TBI ACS/HST has a small field of view, implying 
that the results are quite sensitive to small-scale substructure. It is not un- 
likely, given what we know about the halos of the Milky Way and M31, that 
such a field would be dominated by just one stream such as for example the 
Sagittarius or the Giant stream around Andromeda. It may be misleading 
to correlate the properties of a dynamically young component such as the 
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(outer) stellar halo, with the global propert ies of the paren t galaxy, until a 
more global view of these halos is available ( Ferguso 

The avenue of mapping and resolving the stellar halos of nearby disk 
galaxies is clearly a very worthwhile goal, and one that should be pursued in 
the near future. The advent of Extremely Large Telescopes (ELTs) should 
allow wide-field surveys of the halos of nearby galaxies (also of earlier type), 
perhaps out to the Virgo cluster. A significant addendum would be spec- 
troscopy, both to determine metallicities as well as the dynamics in the out- 
skirts of these objects. In the hierarchical paradigm all galaxies should have 
experienced mergers, and most frequently minor mergers because most of 
the galaxies in the Universe are dwarfs. This implies that their debris should 
end up in a stellar halo, and so this component is expected to be ubiquitous. 
Therefore, unique insights into the assembly histories of galaxies as well as 
on their dynamics can be retrieved by studying their stellar halos. 



8 Epilogue 

The assembly process of a galaxy leaves imprints in the spatial, kinematical, 
age and chemical distribution of its stars. We have argued here that the 
Galactic halo should hold some of the best preserved fossils of the formation 
history of our Galaxy. 

A stellar halo may well be the most common galaxy component and a 
natural outcome of the hierarchical build-up of galactic systems. For example, 
mergers are expected to deposit debris in a spheroidal like configuration. In 
the case of our Galaxy, bo th the halo field population, as well as the globular 
clusters ("Cote ct al."2000|l may h ave such an origi n, as ev idenced by e.g. the 
Sgr dwarf (Bcllaz zini. Ferraro fc Ibata,2003.; Geisler et al.|[2007.1 . Just like for 
stars, it is also possible that some fraction of the globular clusters formed in 
situ, perhaps during early gas-rich m ergers that were presumably extremely 



common in the past ( Kravtsov fc Gn edin 200f|). By studying the dynamical 



properties of the globular cluster population, and linki ng these to those of the 
field stars, such formation scenarios can be tested (e.g. lGasetti-Dinescu et ahl 
120071 ). 

At present it is unclear what fraction of the Galaxy's stellar halo has been 
built via dissipationless merging, and which by the gravitational collapse of 
one or more gas clouds (i.e. via a monolithic-like collapse or during gas rich 
mergers). There is significant evidence of substructure in the form of streams 
in the outer hal o, but this is still at a rather qualitative level (although see 
iBell et al]l2008[ ). It is also very likely that there is currently a bias towards 
detecting the highest surface brightness features, which presumably corre- 
spond to the most recent accretion events of relatively massive satellites. 
Deep wide-field photometric surveys should allow us to quantify the amount 
of substructure in the outer halo in an unbiased manner, especially if ac- 
companied by spectroscopic surveys. Accurate photometry to enable quick 
and easy identification of suitable halo tracers are a must. Planned surveys 
such as VISTAH, the OmcgaCam on the VLT Survey Telescope, SkyMap- 

(http://www.vista.ac.uk, 
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per, PanStarr^ and the Large Synoptic Survey Telescope (LSStQ) will 
significantly contribute to this enterprise. Wide-fie ld imaging of nearby disk 
galaxies as ongoing with HST (Ide Jong et 311120071 ) . and as discussed for the 
ELTs, will allow us to establish the commonality of such merger events. 

The dynamical timescales are much shorter in the inner regions of galax- 
ies, which implies that halo substructures are no longer apparent in physical 
space, but become very prominent in velocity space. For example, a galaxy 
of the size of the Small Magellanic Cloud should have contributed ^ 50—100 
streams /cold moving groups near the Sun if accreted roughly 10 Gyr ago 
(|Helmi fc Whitelll999t) . This implies that full phase-space information is ab- 
solutely vital to recover the ~ 500 streams predicted to be present in the 
inner stellar halo if this had been built up only via "dry" mergers. Current 
datasets are too small, or only have a limited amount of sufficiently accurate 
phase-space information to allow testing such a scenario. Of course this im- 
plies that it is still possible that the halo formed largely during a gas rich 
(merger) phase. 

The definitive understanding of the structure and history of the Galac- 
tic stellar halo is likely to be established only via a "multi-dimensional" ap- 
proach. Knowledge of the stars kinematics complemented with their chemical 
abundances and ages should allow us to address questions such as: 

1. What fraction of the halo has been built via accretion? How docs the 
importance of this process change with distance from the Galactic center? 

2. When was the halo in place? 

3. How many accreted objects have contributed to build up our Galaxy? 
Were these objects major contributors of cold gas, that later produced 
the disk? Is there an evolutionary link between the disk and halo? 

4. What is the relation between the bulge and the stellar halo? The bulge is 
as old as the halo but more metal-rich, and pressumably built partly by 
a dynamical instability which led to the bar. Yet, a fraction of its stars 
may have the same origin as some halo stars. 

5. What were the properties of the accreted objects? What were their masses, 
their star formation and chemical histories? How do they compare to the 
small galaxies we see around us today? 

6. What is the dynamical history of our Galaxy? How did its mass and 
gravitational potential evolve in time? 

The answers to these and many other questions will all be within reach 
in the coming decade, when large datasets with full phase-space and astro- 
physical information become available, i.e. in the era of Gaia. 
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